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Abstract 

We investigate the influence of step edge diffusion (SED) and des- 
orption on Molecular Beam Epitaxy (MBE) using kinetic Monte-Carlo 
simulations of the solid-on-solid model. Based on these investigations 
we propose two strategies to optimize MBE growth. The strategies are 
applicable in different growth regimes: During layer-by-layer growth 
one can exploit the presence of desorption in order to achieve smooth 
surfaces. By additional short high flux pulses of particles one can in- 
crease the growth rate and assist layer-by-layer growth. If, however, 
mounds are formed (non-layer-by-layer growth) the SED can be used 
to control size and shape of the three-dimensional structures. By con- 
trolled reduction of the flux with time we achieve a fast coarsening 
together with smooth step edges. 

PACS: 81.10.Aj; 68. 35. Fx; 68.10.Jy 



1 Introduction 

The growth of high quality compound semiconductors is of great technolog- 
ical importance Despite the longstanding tradition of molecular beam 
epitaxy (MBE), it is still a challenging task to improve the growth of high 
quality thin films and well defined interfaces. In order to optimize MBE 
growth a detailed knowledge of the relation between microscopic processes 
and macroscopic properties is very important. Computer simulations are 
an ideal tool to access this relation between atomistic processes and epitax- 
ial growth. In addition, different and new growth strategies can be easily 
implemented and tested in a fast and cheap way ^. 

In this paper we will investigate the macroscopic effects of two distinct micro- 
scopic mechanisms. The term microscopic refers to processes on the atomic 

^Corresponding author. FAX:+49-931-888-4604; e-mail: schinzer@physik.uni- 
wuerzburg.de 
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scale: e.g. a single diffusion step of an adatom or desorption of an atom. 
These processes are the ingredients to the computer model used in this pa- 
per. This is contrasted to the term macroscopic for effects which are typi- 
cally measurable in experiments: e.g. the overall mass desorption as can be 
monitored by the partial pressure , the form and the distribution of three- 
dimensional structures accessible by scanning tunneling microscopy or 
the growth rate as determined by electron diffraction oscillations [§. The 
computer simulations employed here are ideally suited to bridge the gap be- 
tween such macroscopic effects and their underlying microscopic processes, 
since both scales are accessible. 

Several strategies have been proposed in the literature to optimize MBE 
growth: In particular, layer-by-layer growth is most desirable in order to 
achieve high quality thin films 0, |^. However, quite often a transition to 
non-layer-by-layer is observed where three-dimensional (3D) structures such 
as mounds or pyramids appear. In conventional MBE the time tx until 
this growth mode crosses over to 3D-growth has been shown to vary with 
Fix ~ {D / FY 0, where F stands for the flux and D for the diffusion 
constant of adatoms. Without desorption, Ehrlich-Schwoebel barriers and 
step edge diffusion (SED) 5 = 2/3 has been observed for epitaxial growth 
For metals, several methods have been proposed and tested to achieve and 
maintain layer-by-layer growth. For instance it has been shown that pulsing 
the deposition rate or pulsing the temperature leads to a prolonged layer-by- 
layer regime pf. Recently, it has been proposed that pulsed glancing-angle 
sputtering can even lead to "layer-by-layer growth forever" |]^. All these 
concepts can so far be understood in terms of a typical diffusion length or 
an enhanced interlayer diffusion at step edges. 

In this paper we will propose strategies which exploit other specific micro- 
scopic processes, namely desorption [1^ and SED [11, 1^. As far as we know. 



no attempt has been made to exploit these processes in order to achieve 
improved growth. Some preliminary results of our investigation have been 
published in , and in this paper we describe the investigation in full detail. 

In Sec. we introduce the sohd-on-solid (SOS) model and the microscopic 
processes. In our computer experiments we first investigated the tempera- 
ture dependence of the overall growth rate in the layer-by-layer regime (Sec. 

We are able to correlate this macroscopic property to the microscopic 
dynamics of the computer model. This allows us to propose a new strat- 
egy for layer-by-layer growth. If, however, the growth of three-dimensional 
structures occurs another strategy is applicable. Using a simplified model of 
growth we have recently shown that SED plays a crucial role in this regime 
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JL^ , 0. These findings allow us to propose an optimized way for the growth 
of SD-structures in Sec. ^ Concluding remarks concerning the experimental 
realisation and a summary will be given in Sec. |^. 



2 Computational model 

Lattice models with the solid-on-solid restriction (SOS) have been proved 
to be a useful tool to study surface morphology ||16|, The model has a 
long history for the study of the surface roughness transition |TB[. Gilmer 
and Bennema were the first (to our knowledge) who included surface diffusion 
]. Since then it has been intensively used to study epitaxial growth 



Here we use its most simple form where only one kind of particles and a 
simple cubic lattice is considered. The particles represent single atoms when 
a comparison with a simple cubic metal is made. However, even compound 
semiconductors can be modeled, as long as kinetic features are investigated 
only. E.g. in Ref. the RHEED (reflection high energy electron diffrac- 
tion) oscillations of GaAs(OOl) during growth have been quantitatively re- 
produced. 



In our simulations we use the Maksym-algorithm of Ref. . At each time 



step a Monte-Carlo move is carried out. The way how the event is selected 
makes it superior to conventional Monte-Carlo techniques (the algorithm 
uses partly a binary search in the array of possible events). We have used a 
system of 300 times 300 lattice sites, if not stated otherwise. 

Besides the SOS-restriction further simplifications are due to the particular 
choice of possible events labeled i and the parametrisation of the correspond- 
ing rates Fj. We allow jumps to the four nearest neighbor sites (diffusion on 
a flat surface, attachment and detachment from steps, . . . ) and desorption. 
The rates do only depend on the four nearest neighbor sites as will be de- 
scribed below. We describe all these processes as Arrhenius-activated 

r, = ..exp(-^), (1) 

as is predicted by several theories p4| . 



One quite often assumes vibration frequencies Ui of the order of Debye- 
frequencies, i.e. 10^^ - lO^^s"^. Indeed, in sublimation experiments of 
CdTe(OOl) lO^'^s"^ has been observed vibration frequencies for dif- 

fusion are often of the order of lO^^s"^ (measurements for metals ||26|, cal- 
culations for GaAs(OOl) or simulations and calculations for Si(OOl) 
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29], 0). Hence, it is reasonable to assume that the diffusion as well as the 
desorption rates of our model share one common prefactor Ui = = lO^^s^^ 
which allows to keep the number of parameters small. 

The activation energy for the different microscopic processes are parameter- 
ized as follows: a diffusion jump of a free adatom has to overcome a barrier 
Eb, each next in-plane neighbour adds an energy Ei\i. The rate of diffu- 
sion jumps which keep the height of the particle unchanged thus becomes 
uo exp {—{Eb + nE^) /k.^T), where n represents the number of next in-plane 
neighbors. Note, that the overall rate for diffusion on a fiat surface is four 
times this jump rate due to the four possible directions. Hence, the diffusion 
constant becomes D = i'Qexp{—EB/kBT) [^. Since we measure all length 
scales in units of the lattice constant a we have neglected the term in 
D. At step edges an additional Ehrlich-Schwoebel barrier Es is considered 



32| , p3|| . However, this barrier is not added for particles on top of an elon- 
gated islands of only one lattice constant width The desorption barrier 
is Eb). Again, each next in-plane neighbour contributes Ejsf. 

The deposition of particles occurs with a rate F measured in monolayers 
per second (ML/s). During deposition we consider another process which 
is not Arrhenius-activated. After a deposition site is chosen randomly we 
allow the particle to relax to a lower neighboring site. Here, we consider only 
relaxation to nearest-neighbor sites. Such transient diffusion or downward 
funneling has been observed in molecular dynamics of simple Lennard- Jones 



systems p5|, B6| and has been related to the reentrant layer-by-layer growth 



at very low temperatures fS^]. In addition it has been shown to play a crucial 



role for slope selection in mound morphology [|I^ . 

We will concentrate on one set of parameters, namely Eb = 0.9eV, E]\f = 
0.25eV, Es = O.leV, and Eb = l.leV. This particular choice of parame- 
ters reproduces some features of CdTe(OOl) during sublimation [^, ^ and 
annealing. However, we would like to stress that the findings of our present 
work are of more general relevance, independent of the specific choice of the 
energetic parameters. 



3 Reevaporation during layer— by— layer growth 
and the Flush Technique 

For clarity we will distinguish several processes of desorption. The term 
desorption will be explicitly used to describe the atomistic process: the des- 
orption of a single atom. Sublimation is reserved to describe the evaporation 
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Figure 1: Reevaporation rate during growth with F = 1 ML/s (O), F = 
4 ML/s (O), and subhmation rate xlO at F = ML/s (A). In addition 
we show the reevaporation rate using the proposed flush-technique (□) with 
a mean flux of 1 ML/s consisting of a constant flux of 0.77 ML/s plus an 
additional pulse of 0.23 ML during 0.003 s at the beginning of each second. 



of a surface when left in (perfect) vacuum. Reevaporation or more precisely 
reevaporation during growth will describe the overall desorption rate mostly 
due to the desorption of freshly deposited particles during growth. 

For two of the most important compound semiconductors a decrease of the 
MBE growth rate with increasing temperature was observed (CdTe(OOl) 
P, H, H and GaAs(OOl) [|l|, H). For CdTe(OOl) the reevaporation rate 
was found to follow an Arrhenius-rate with considerably lower values of the 
activation energy of 0.14 - 0.30 eV compared to sublimation (1.55 eV resp. 
1.9 eV 1^). A tempting explanation is to ascribe this low energy to the 
existence of a physisorbed precursor |2^. However, studies of the sublima- 
tion with computer simulation as well as experiments for CdTe(OOl) ||39|1 
showed a strong influence of the morphology. In Ref. PB[ we already con- 
cluded that in MBE, one should expect desorption rates other than those 
measured by sublimation. Independently Pimpinelli and Peyla also showed 
that a physisorbed precursor is not necessary to explain the observed low 
energies using kinetic Monte-Carlo simulations as well as simple scaling ar- 
guments [||, I]]. 



In flg. H, the diamonds represent the reevaporation which we derived from 
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Figure 2: Variation of the reevaporation rate during growth as a function 
of the desorption barrier Ed (upper) and diffusion barrier Eb (lower curve). 
The slope gives the contribution to the effective energy. The simulations 
were carried out at T = 560 K with F = 4 ML/s. 



the difference between the applied flux {F = 1 ML/s) and the measured 
growth rate {i.e. the reached height /simulated time). The data points for 
F = 4 ML/s (O) show that the effective energy is independent of the applied 
flux. The triangles correspond to the sublimation i.e. the evaporation rate 
without application of an external flux {F = ML / s) . Both processes are 
found to be Arrhenius-activated, however, with strikingly different effective 
energies. The reevaporation rate during growth corresponds to an activation 
energy of approximately 0.90 eV which is even lower than the microscopic 
desorption energy of E^ = 1.1 eV. At high temperatures the reevaporation 
rate saturates and equals the flux of impinging particles. On the contrary, 



the sublimation energy of approximately 1.73 eV is considerably higher |45 



The relation of the sublimation energy to the microscopic parameters was 
shown to be approximated by E^nh ~ O.QIEb + 0.35Ed + 2.85En + 0.44^^ 
P5| . To derive this relation we varied all microscopic energy parameters 



independently. Applying the same microscopic analysis to the reevaporation 
during growth of this model we obtain 



E ^ 



-O.SIEb + 0.94Ed + 0.51En - 0.03^^. 



(2) 



As an example for this microscopic analysis, fig. ^ shows the measured in- 
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Figure 3: Schematic variation of the flux using the flush-technique. 

fluence of the diffusion barrier Eb and of the desorption barrier Ej:, to the 
reevaporation rate. We apphed a flux of F = 4 ML/s. Since the measured 
reevaporation rate is much lower (less than 1 ML/s) we can be sure to have 
no saturation effects. Note the opposite sign of the two contributions. The 
slope measures the prefactor in the above expression of E-^^, {eq. |^). 

We want to mention that this result does not agree with the scaling relation 



obtained by Pimpinelli and Peyla [Q. However, at lower temperatures (not 



shown) we observe a cross-over to their result with a critical nucleus size of 
i* = 1. The crossover itself can be seen at the data point for F = 4 ML/s at 
500 K which lies above the value of about 0.2 ML/s which is extrapolated at 
500 K from data points at higher temperatures. A detailed investigation of 
the validity-regime of our result and the relation to the results obtained by 
Pimpinelli and Peyla will be postponed to a future work. 

Besides the different weightings in E^.^ and E^nh the striking difference (at 
high as well as low temperatures) is the negative contribution of the diffu- 
sion barrier Eb to E^^- This result seems to be of general validity |^ and 
can be explained in the following way: Even though the island distance is 
influenced by Eb, the dominant effect of higher diffusion barriers seems to be 
the reduction of the diffusion length of the adatoms. Consequently, particles 
have a higher probability for desorption before they stick to an island. 

This result suggests a strategy to obtain high-quality (layer-by-layer) growth 
together with high growth rates. Short flushes of particles at the beginning 
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Figure 4: Comparison of the surface width for conventional growth (A) 
and the flush-technique as described in the text. We consider the flush- 
technique with (B) and without desorption (C). The inset shows the surface 
width oscillations during the deposition of the first 20 ML. For clarity of 
presentation we have shifted the upper (lower) curve about +0.2 (-0.2) inside 
the inset. These simulations were carried out on a 512 x 512 lattice at 560 K. 
The results are averaged over six independent simulations runs in each case. 

of each monolayer would result in a great density of islands. Afterwards with 
a low flux the particles probably hit islands to stick to which will result in 
a low overall reevaporation rate. The proposed procedure (flush-mode) is 
drawn schematically in fig. ^ 

Figure |l| shows that the reevaporation rate indeed reduces by a factor of 
about two when applying this strategy. The mean flux was 1 ML/s as for 
the conventional growth simulations. The profile of the flux was composed 
as follows: In intervalls of one second we deposited a total amount of 0.23 
ML within 0.003 s (see fig. ^). Afterwards a constant flux of 0.77 ML/s was 
applied. According to the decrease of evaporation the growth rate increases. 
The gain is highest at high temperatures (at 620 K the growth rate is doubled) 
since there the evaporation rate becomes comparable to the applied flux. 

In addition to higher growth rates, layer-by-layer growth is assisted by 
the flush-mode. In fig. ^ we compare three different techniques/models of 
growth: 

A) conventional growth with F = 1 ML/s and allowed desorption, 
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B) flush-mode with Fconst = 1 ML/s and additional 0.30 ML in 0.003 s 
each second and allowed desorption, 



C) flush-mode without desorption [Ed — oo), -Fponst ~ 

0.77 ML/s and 

additional 0.23 ML during 0.003 s each second. 

The different fluxes in B and C are chosen in order to achieve the synchro- 
nization of the pulses with layer-completion. Due to the possible desorption 
in B, however, synchronization can be achieved only approximately in this 
case. 

We investigate these different methods by comparing the surface width 



Perfect layer-by-layer growth would lead to oscillations between zero and 
0.5 (coverage of half a monolayer). Higher values of w indicate a broader 
distribution of the heights. 

After the deposition of 60 ML using the different techniques the surface 
widths become considerably different (see fig. §). The flush-mode without 
desorption (C) is even farer away from perfect layer-by-layer growth com- 
pared to conventional growth (A). The flush-mode in the presence of des- 
orption (B) is superior to both, (A) and (C) in the long run and keeps the 
surface smooth. Looking at the deposition of the first 20 ML this seems to be 
surprising. The oscillations of w with technique (B) are disrupted due to an 
obvious asynchronization. With (C) the synchronization is perfect leading 
to very strong and regular oscillations. Using the conventional technique (A) 
the oscillations are damped and much less pronounced. 

To summarize, the usage of the proposed flush-technique is useful to improve 
the growth rate (we achieved a factor two at high T) and to assist layer- 
by-layer growth. Hereby, the desorption of adatoms is crucial to achieve 
optimized growth: without desorption the flush-mode is worse compared 
to conventional growth even though strong oscillations are induced. The 
reevaporation has such an impact because it is height selective, i. e. adatoms 
on top of existing islands desorb easily whereas adatoms beneath islands 
preferentially are incorporated in the crystal. Clearly, this height selective 
behaviour is achieved only when a positive Ehrlich-Schwoebel barrier hinders 
the particles to be incorporated at step edges from above. 

We would like to point out that the usage of a chopped flux has been proposed 
and investigated by Rosenfeld et al. in the framework of the concept of 




(3) 
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the two mobilities. However, our findings show that only in the presence 
of desorption the occurrence of oscillations is indeed coupled to a reduction 
of the surface roughness. In ref. the effect of a chopped flux on island 
distances was investigated. These findings would allow to optimize even 
further in that one calculates the minimal flux intensity and the time of the 
flush needed in order to achieve an increased island density. Here, we have 
chosen a safe high flux without explicit use of the results of ffl . 



4 Optimizing the structure of mounds in 3D— 
growth 

Quite generally, layer-by-layer growth as well as step flow is not attain- 
able forever This can be due to e.g. Ehrlich-Schwoebel barriers 



49 1 which is typically positive |50]. This favors new nucleation events on 
top of existing islands which leads to 3D-growth sooner or later. In order 
to optimize MBE growth it is thus also interesting to study the growth of 
SD-structures by computer simulations. 

We will start with a brief summary of our findings for 3D-growth on the 



basis of a simplified model of epitaxial growth [|1^, |T5[. This will enable us 
to introduce the basic concepts. After that, we will show how these results 
can be used in order to improve 3D-growth (to be specified below). We will 
test this new strategy with computer simulations of the SOS-model of Sec. 



The most important simplification we introduced in ref. |T^ was an effec- 
tive description of diffusion and nucleation. Rather than to simulate the 
simultaneous motion of many adatoms we concentrated on the simulation 
of individual particles which is a usual technique for simple growth models 
Parameters to the model are the diffusion length and in a similar 



way SED is considered. Even though a similar SED has been introduced in 
ref. [^, there, in difference to the present work no search for kink sites was 
implemented. 

In MBE the typical length of the step edge diffusion process depends on tem- 
perature and flux of arriving particles . On a one-dimensional substrate 



the theory of island nucleation predicts a typical distance between nucleation 
centers of the form 

isED ^ h (4) 
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where d is the diffusion constant and / the ffux of arriving particles [Q]. If 
we apply this theory to the lateral or in-plane growth of a pyramid (concen- 
trating on a slice of one ML thickness) the flux / can be identified with the 
reduced flux per unit length of the step edge / = Fix where It stands for 
the terrace width. Within this context d becomes the diffusion constant for 
diffusion along the step edge. The scaling relation (|D for ^sed was obtained 
under the restriction that two atoms (i.e. i* + 1 = 2) already form a stable 
nucleus and no desorption occurs. We note that for greater values of i* the 
correct theoretical result has been derived recently ||5^. However, for the 



model as described in Sec. ^ the assumption of = 1 is reasonable. 

When £sED is of the order of the modeled system size (strong SED) the 
growth is characterized by the formation of square based pyramids with a 
well defined slope. The step edges are oriented along the lattice coordinates 



and the surface width was found in Ref. |T3 to grow with a power law 



w och'^ with (3 ^ 0.45 (5) 

where P is called the growth exponent |T^. Typically, one expresses the 
scaling behaviour in terms of the elapsed time. However, in the context of 
this paper it is advantageous to use the mean height h instead (as will become 
clear soon). The typical distance between the pyramids (the correlation 
length) was found to be proportional to 

^ oc h^/' and z = a/P^ 2.3 (6) 

in accordance to the occurrence of slope selection. More formally this means 
that the ratio of the typical length scales, w and ^, remains constant, and 
hence a = 1. 

The relatively high growth exponent of 0.45 reflects the fact that the coars- 
ening process is SED-assisted [|T^]. Due to the strong SED material is moved 
efficiently towards regions with high densities of kink sites, i. e. towards the 
contact points of pyramids or mounds. For lower values of ^sed the coarsen- 



ing process is purely noise assisted |5^. Hence, the structures are merging 
more slowly. 

If the size of the pyramids exceeds £sed the pjnramids loose their perfect 
shape. The structures become round, and step edges will be fringed []T3|. It 
is clear that due to the coarsening, conventional MBE growth is bound to 
drive itself into this state. 

Now, we turn to the investigation how the latter stage in MBE growth can 
be prevented. The main idea is that in order to prevent the occurrence of 
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Figure 5: Flux variation used in the simulation of flux-adaption of fig. 
The inset shows the reached height in dependence of the time. 

rough step edges one has to require that always £sed ~ C In the follow- 
ing we demonstrate how to fulfill this condition by varying the flux F of 
arriving particles. Equally well, one could adapt the growth temperature. 
However, for this a detailed knowledge of the activation energy of d and the 
temperature dependence of It IS^ is necessary which is often not available. 



Equating the expressions and (H) we obtain the height dependence of the 
flux 

F{h) = ch-^/' (7) 

where c is an adequate constant. To reformulate this relation in terms of 
the time we use dh/dt = F and solve the resulting differential equation, 
obtaining 

h{t) oct-'/^^+'\ (8) 

Reinserting this result into (|^ we obtain that the flux should be varied 
according to 

F{t) oc t-4/(4+^) ^ (9) 

where we inserted z = 2.3 according to SED-assisted coarsening [p^ . 

We applied this strategy to the growth of the SOS-model of Sec. at 560 
K. Clearly, SED is not a process which is explicitly considered in this model. 
Typically, the atoms with only a single bond to a step edge will detach and 
diffuse on the terrace. However, the net result will be the same: single 
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Figure 6: Comparison of the surface width under optimized growth con- 
ditions (A) and without an adaption of the flux F = 1 ML/s (O). Under 
optimized growth conditions the flux was initially set to 2 ML/s and after 
10 s adapted according F = 2/(t/10s)°-^^ML/s (see fig. 1). The sohd lines 
correspond to growth exponents (5 = 1/3 and f5 = 1/2. The arrows mark the 
positions of the snapshots of fig. |^ 



bonded atoms will be moved to places with higher coordination (kink sites). 
To prevent the inference of reevaporation we suppressed this process ||58|| . 
However, we checked that even with desorption, the strategy is still applicable 
and useful. The flux was chosen as shown in fig. ^ We started with a constant 
flux of 2 ML/s. After the growth of twenty monolayers we adapted the flux 
with time t according to F = Fo/(t/10s)°'^^. 

In fig. 1^ we compare the resulting evolution of the surface width w with 
a simulation with a constant fiux F = 1 ML/s. With the new strategy 
we obtain a higher growth exponent of /3 ~ 1/2 compared to conventional 
growth with j3 ~ 1/3. These exponents fit well to (3 = 0.45 for strong 
SED |l^ and l3 = 0.33 for intermediate values of ^sed HH- The result 



that w grows fast and is described by a high growth exponent under these 
optimized growth conditions should not however be confused with the notion 
of a fast roughening, self-affine surface. It just means that the structures are 
merging fast and the mounds are getting high and wide. This can be seen 
directly in fig. |^ After the deposition of 300 ML under constant fiux the 
structures are small whereas under optimized growth conditions the resulting 
structures are larger. Note that because of the higher initial fiux of 2 ML/s 



13 



Figure 7: Snapshots of surface morphology of 300 x 300 lattices where 
about 300 ML were deposited with the two different methods as in fig. ^ 
(left: flux-adaption, right: conventional growth). 

the island density was much higher in the beginning under optimized growth 
conditions. Nevertheless the SED-assisted coarsening leads to a considerably 
fewer number of mounds (approximately 10 which should be compared to 20 
with the conventional growth). 

It is clear that in order to obtain larger structures in conventional MBE one 
just has to grow for longer times. The step edges will become smooth due 
to the equilibration after the growth has been stopped. However, during 
growth the step edges do not remain smooth as in our optimized growth 
mode. Hence, a larger probability for the creation of vacancies or other 
crystal faults will be present. After growth stops these faults can probably 
be only partially eliminated in the non-optimized growth. We also mention 
that in the end another processes will of course become important too. In 
the limit of t — oo no net growth will be achieved in the optimized growth 
and the equilibration of the surface will be dominant. 

5 Conclusion 

We have investigated the effect of the microscopic dynamics on experimen- 
tally accessible macroscopic effects in MBE growth. Based on simulations of 
the solid-on-solid model we proposed two optimized growth-strategies. 

Comparing the layer-by-layer growth with sublimation we understood how 
the desorption of single adatoms comes into play during growth. During 
growth, freely diffusing adatoms are created by the external flux. During 
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sublimation, however, such adatoms must first be created, e.g. through de- 
tachment from steps. This difference manifests itself in the different contri- 
bution of the microscopic activation energies to the effective energies. The 
diffusion barrier Eb increases the effective energy of sublimation whereas it 
decreases the activation energy of the desorption rate during growth. Since 
the macroscopic desorption rate is influenced by the typical lifetime of single 
atoms we are able to intervene in the growth process. We showed that a flush- 
mode is able to prolong the layer-by-layer growth regime and to reduce the 
desorption rate: applying short pulses of particles we create a high density 
of islands. Afterwards with a low flux one completes the monolayer. At least 
for our particular simulations the desorption was crucial to obtain improved 
growth. Even though the flush-mode always induces strong oscillations, only 
in combination with desorption it leads to an improved growth. This can be 
explained with the height-selective behaviour of desorption, i.e. desorption 
occurs preferentially on top of islands as long as a positive Ehrlich-Schwoebel 
barrier is present. 

In experiments one should be able to produce such short flushes using a 
chopper or pulsed laser deposition in conjunction with conventional MBE. 
The method should be very useful in order to grow planar coherent thin fllms, 
e.g. for application in quantum well structures. In addition, such experiments 
would allow to decide wether desorption occurs out of a physisorbed precursor 



which has been debated in the literature 25, 421. If so, one should not obtain 



an improved growth rate using the flush-mode. 

However, since layer-by-layer growth is unstable, the optimization of 3D- 
growth might be useful too. Based on recent results concerning the step 
edge diffusion we proposed to vary the flux of arriving particles in order to 
maintain smooth edges already during growth. Reducing the flux according 
to F oc we were able to recover the high growth exponent of /5 ~ 0.45 



measured on the simplifled MBE model |T^. The fast coarsening process 
(since SED assisted) yields structures which soon become very large com- 
pared to those of conventional MBE. Irrespective of the desired size of the 
structures they can be produced under the same (strong SED) conditions 
which is accomplished by variation of F. Otherwise the MBE growth would 
drive itself in the regime where ^sed is less then the typical extension of the 
structures. Thus, our method opens new possibilities for the controlled cre- 
ation of these selforganized nanostructures by MBE. In addition, this strategy 
should reduce the probability for the creation of vacancies since during the 
conventional growth the rough edges would be overgrown later. However, this 
is speculative and cannot be verifled in the framework of the solid-on-solid 
model. 
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Typically, rather low fluxes are used in order to improve the quality of the 
grown structures. However, our result suggests that it is not disadvantageous 
to apply higher fluxes in the beginning. In the end, when the resulting 
structures arc rather large it becomes important to reduce the flux in order 
to adapt it to the smoothening range of the step edge diffusion. 

*** 
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